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Modeling the chemical evolution of to Centauri using 
three-dimensional hydrodynamical simulations* 
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ABSTRACT 

We present a hydrodynamical and chemical model for the globular cluster us Cen, 
under the assumption that it is the remnant of an ancient dwarf spheroidal galaxy 
(dSph), the bulk of which was disrupted and accreted by our Galaxy ~10 Gyr ago. We 
highlight the very different roles played by Type II and Type la supernovae (SNe) in 
the chemical enrichment of the inner regions of the putative parent dSph. While the 
SNe II pollute the interstellar medium rather uniformly, the SNe la ejecta may remain 
confined inside dense pockets of gas as long as succesive SNe II explosions spread them 
out. Stars forming in such pockets have lower a-to-iron ratios than the stars forming 
elsewhere. Owing to the inhomogeneous pollution by SNe la, the metal distribution of 
the stars in the central region differs substantially from that of the main population of 
the dwarf galaxy, and resembles that observed in us Cen. This inhomogeneous mixing 
is also responsible for a radial segregation of iron- rich stars with depleted [a/Fe] ratios, 
as observed in some dSphs. Assuming a star formation history of ~1.5 Gyr, our model 
succeeds in reproducing both the iron and calcium distributions observed in us Cen 
and the main features observed in the empirical a/Fe versus Fe/H plane. Finally, 
our model reproduces the overall spread of the color-magnitude diagram, but fails in 
reproducing the morphology of the SGB-a and the double morphology of the main 
sequence. However, the inhomogeneous pollution reduces (but does not eliminate) 
the need for a significantly enhanced helium abundance to explain the anomalous 
position of the blue main sequence. Further models taking into account the dynamical 
interaction of the parent dwarf galaxy with the Milky Way and the effect of AGB 
pollution will be required. 
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1 INTRODUCTION 

The stellar system us Cen (NGC 5139) is unique amongst 
Galactic star clusters in terms of its structure, kinemat- 
ics, and stellar content. It is the only kn own globular clus- 
ter (GC) showing a clear [Fe/H] spread (jNorris et all 1 19961. 
and references therein). Recent photometric surveys have 
revealed the presence of multiple sequences in its color- 
magnitude diagr am (CMP), indica t ing a complex star for- 
mation history (ICannon fc Stobiel 1 1973! : iRev et all |2004 
ISollima et al.ll2005ah . On the basis of the analysis of the red 
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giant branch (RGB) morphology, three different stellar com- 
ponents have been identified: a dominant metal-poor popu- 
lation ([Fe/H]~ —1-6), an intermediate population spanning 
the metallicity range —1.3 <[Fe /H]< —1.0, and a me tal rich 
component with [Fe/H]~ -0.6 l|Pancino et al.ll2000h . 

Peculiarities have also been found along the Main Se- 
quence (MS) of the cluster where an additional blue MS 
(bMS, comprising ~ 25% of the cluster's MS stars) run- 
ning p arallel to the dominant sequence ( rMS) has been re- 
solved ([Anderson 2002; Bcdin et al. 2004). In contrast with 
that predicted by stellar models with canonical chemical 
abundances, bMS stars show a metallicity a fact or of two 
higher than that of the rMS l|Piotto et all 120051 ). A large 
helium overabundance (AY ~ 0.15) in the bMS sequence 
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has be en proposed to explain its anomalous position in th e 
CMD l|Norrisl 12004 iPiotto et ai]|2005l ; ISollima et al.ll2007l ). 
However, such a large helium abundance spread poses se- 
rious problems for the overall interpretation of the chem- 
ical enrichment history of this stellar system. In fact, no 
chemical enrichment mechanism is able to produce the huge 
amount of helium required to reproduce the observed MS 
morphology wit hout dramatica l ly imp a cting upon the metal 
abundance (cf. iKarakas et al.l [2006a: Soll ima et al.l 120071 : 
iRomano et al.ll2007f ). ~ 

In terestingly, IPancino et ail (|2002h and lOrigl ia et al.l 
i|2003l ) found that while the metal-poor and intermediate- 
metallicity stellar populations of uj Cen have the expected 
a-element ove rabundance observed in halo and GC stars 
([a/Fe]) ~ 0.3 (lEdvardsson et al.ll 19931 ). the most metal-rich 
population ([Fe/H]~ —0.6) shows a significantly lower a- 
enhancement (([a/Fe])=0.1). The chemical composition of 
the former populations requires that at least part of the 
gas released in the interstellar medium (ISM) by SNe II 
must have been retained by uj Cen over a relatively short 
(<1 Gyr) time interval. However, the presence of stars with 
lower values of [a/Fe] at larger values of [Fe/H] indicates 
SNe Ia have also likely contributed to the chemical evolu- 
tion of uj Cen. Since SNe Ia explosions occur over a timescale 
longer than that of SNe II, the ratio [a/Fe] of the polluted 
ISM tends to decrease with time as [Fe/H] increases. Al- 
though the chemical properties of the stellar populations 
of uj Cen indicate a prolonged star formation, there is no 
consensus as to its duration. An age spread of 3-5 Gyr has 
been claime d to explain the overa ll trend of a-elements with 
metallicity l|Romano et al-lfeoOTh . A somewh at shorter (< 
1.5 Gyr ) age difference has been suggested bv lSollima et al.l 
(|2005bl ). comparing theoretical isochrones with the location 
of a sample of sub-giant stars with known metallicity. 



In any case, the body of evidence collected so far leads to 
the hypothesis that uj Cen was formerly a larger stellar sys- 



tem, possibly a dwarf galaxv (|Hughes & Wallersteinl 


2000; 


Dinescu et al.lll999; Maiewski et al.ll2000l;ISmith et all 


2000; 


Gnedin et all 20021; Bekki & Norrisl 20061: Romano et al. 


20071). that lost most of its stars and Eas in the interaction 



with the Milky Way. Under this hypothesis, self-consistent 
dynamical m odels, and coupled N - body + hydrodynamical 
simulations dCarraro fc Lial |2000| ; iBekki fc Freemanl 120031 ; 
iTsuchiva et al.l I2004T ). succeed in reproducing the main ob- 
servable characteristics of this system, assuming a total mass 
for the initial object in the range 10 8 -10 9 M . Tidal inter- 
action with the Galaxy, as well as ram pressure stripping 
by its halo gas, have also been invoked to explain the stel- 
lar structure and kinematics of the local d warf spheroidal 
galax i es (dSphs), and th eir lack of gas (e.g. Ivan den Berghl 
Il994l ; iMaver et all 120061 ). In this scenario uj Cen could re- 
tain the ejecta of previous generations of stars and thereby 
self-enrich throughout its star formation history (SFH). 

We propose a model for uj Cen based upon 3D hy- 
drodynami cal simulations of an evolving dwarf spheroidal 
galaxy (see iMarcolini et al.ll2006l ). In Section 2, we describe 
the evolution o f uj C en in the general frame proposed by 
IMarcolini et alj (|2006l ). The general results and the possibil- 
ity of an extended SFH are discussed in Section 3. Section 4 
is devoted to the simulation of the systems' CMDs, while 
our conclusions are drawn in Section 5. 



2 THE MODEL 



2.1 Qualitative framework 

To fully appreciate the chemical enrichment history of uj Cen 
it is necessary to understand both its SFH and the associated 
role of SNe feedback on its dynamical and chemical evolu- 
tion. IMarcolini et alj (|2006l ) provide a useful framework in 
which to pursue work of this nature, in that their three- 
dimensional hydrodynamical simulations of isolated dSphs 
are directly comparable to the objects from which uj Cen has 
been proposed to have originated. In this study, the authors 
assumed a prolonged intermittent SFH, tracking the roles of 
both SNe Ia and SNe II in the chemical enrichment history 
of the system. These models, although intended to give a 
general picture of dSphs, were tailored to the Draco dwarf 
galaxy. 

The total energy released by the SNe II explosions is 
much larger than the binding energy of the gas. However, 
efficient radiative losses enables the galaxy to retain most 
of its gas, providing the necessary fuel source for the afore- 
mentioned prolonged SFH. 

The SNe II are more concentrated toward the galactic 
centre, where their remnants overlap forming a single cavity 
composed of a network of tunnels filled by hot rarefied gas. 
The bulk of the ISM is pushed outward to the edge of the 
cavity. Once the SNe II cease exploding (~30 Myr after each 
star burst episode), the global cavity collapses and the ISM 
flows back down the potential well, giving rise to a new 
starburst (see Fig. 3 in Marcolini et al. 2006) . In conclusion, 
the gas content in the centre has an oscillating temporal 
profile whose period is given by the time interval between 
two successive starbursts. 

Given their lower rate, SNe Ia do not significantly af- 
fect the general hydrodynamical behaviour of the ISM, but 
their role is relevant for the chemical evolution of the stars. 
Because of the longer timescales over which they contribute, 
the SNe Ia progenitors produced in previous starbursts con- 
tinue to explode during the quiescent periods, even after the 
gas has settled back to the central region. The higher am- 
bient gas density (together with the lower explosion rate) 
results in SNe Ia remnants being isolated from one another 
during this phase (see Fig. 14 in Marcolini et al. 2006). As a 
consequence of this inhomogeneous pollution, stars forming 
in the regions occupied by SN Ia remnants (hereafter re- 
ferred to as "SNe Ia pockets") have lower [a/Fe] ratios and 
higher [Fe/H] ratios than those formed elsewhere. This effect 
is particularly important for the chemical evolution of the 
central region, where the SN Ia rate is greater and where 
SN Ia pockets formed in the outer zones are returned to 
the galactic centre, with the gas flow during the re-collapse 
phases. The cycle of expansion and re-collapse experienced 
by the galactic gas tends to homogenise the ISM and SNe 
ejecta rather rapidly (~10 8 yr). 

These authors also found that only ~ 20% of the metals 
ejected by SNe were present in the region where the stars 
were assumed to form, despite the absence of a galactic wind. 
The missing metals were pushed to larger distances by the 
continuous action of the SN explosions and did not enrich 
the forming stars (thus mimicking an "outflow"). 
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Figure 1. [Fe/H] distribution function of the long lived stars for the whole "entire dSph volume" (left panel; r* jC = 130 and r*,t = 650 
pc), and two more centrally located regions called "central region" (middle panel, r* t t = 90 pc) and "nuclear region" (right panel; 
r* t = 30). The dashed lines represent stars with [o/Fe]^ 0.2. See text for more details. 



2.2 SNe chemical abundances 

The numerical method used to treat the SNe explosions 
is described in iMarcolini et alj (|2006l ). Here we describe 
briefly the nucleosynthesis prescription assumed in this pa- 
per. Nucleo synthesis products of SNe II and SNe la are 
taken from llwamoto et alj (1 19991 . and references therein). 
Their Table 3 summarises the mean SN II yields averaged 
over the 10 — 50 Mq precursors, and the yields for differ- 
ent models of SN la (we will refer to their W7 model). We 
calculate the a-element enrichment by summing the differ- 
ent contributions of O, Mg, Ne, Si, S, Ca, Ti and Ar. We 
find that a slightly enhanced Ca production (50%) over the 
mean SNe II value helps to better reconcile model predic- 
tions with observations. Note that this increase lies within 
the uncertainties associated with the different SNe models. 
These yields depen d slightly on metallicity, unless metals are 
comp letely absent l|Wooslev fc Weaverll 19951 : iLimongi et al.l 
2000); we do not consider this dependence because it does 
not affect the [a/Fe] feature s of the stars with [Fe/H] > —3 
jGoswami fc Prantzos|[2"000l ). The adopted yields are listed 
in Table 1, t ogether with the adopte d solar abundances 
(adopted from lGrevesse fc Sauval|[l998h . 



2.3 Color- magnitude diagrams 

As a first step, in order to simulate th e morphology of the 
CM P we use the evolutionar y tracks of ICassisT et ail (|2004T ) 
and lPietrinferni et ail (120061 ). calculated for two different en- 
enhancement levels ([a/Fe]= 0.0 and [a/Fe]=+0.4). As we 
will show in Sec. [3] a number of stars in our models show 
[a/Fe] values that lie outside these limits, especially at sub- 
solar values of [a/Fe]. In order to take these stars into ac- 
count, we extrapolated for each metallicity Z the evolution- 
ary tracks down to values of [a/Fe] = — 0.4. 

For each star, our code stores its iron abundance [Fe/H], 
its metallicity Z, its a-enhancement [a/Fe], and the epoch 
of formation. Then, for each star we randomly extracted 
100 masses from a Salpeter (1955) IMF in the mass range 
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M Fe 
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M Si 
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M 0a 


8.80 x 10~ 3 


1.19 x 10~ 2 
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2.34 


0.43 


1.12 xlO" 2 


[a/Fe] 


0.45 


-1.20 






a From Iwamoto et al. 1999 

b Averaged over the progenitor mass range 10-50 Mq 
c Model W7 

d Solar abundances from Grevcssc & Sauva 3 (119981 s ). 



0.5 Mq <M< 0.9 Mq and placed them in the theoretical 
log(L/L Q ), logT c ff plane by means of suitable interpolations 
of the adopted evolutionary tracks. Luminosities and effec- 
tive temperatures were transformed into the desired photo- 
metric system by interpolation within appropriate tables for 
photometric conversions (for a detailed description see Cas- 
sisi et al. 2004). Following these prescriptions we simulated 
a CMD of 100,000 stars in the Johnson- Cousin and in the 
ACS VEGAMAG photometric systems. 

In section [4] we w ill compar our synthe t ic CMD with 
the ob serve d CMDs bvlSollima et ail (l2005al ). lFerraro et all 
\2004 ) and ISollima et ail (|200^ .~Tb convert the absolute 
magnitudes into the apparent ones, we need to assume a 
distance modulus and a reddening correction. In the fol- 
lowing, we adopt (m - M ) = 13.70 (|Bellazzini et alj|200j; 
Del Principe et al][200rj) . the reddening E(B-V)=0.11 (|Lubl 
20021 ). the extinction coefficients A B = 4.1 E(B - V) 
Av = 3.1 E(B-V), A R = 2.32 E(B-V) l | Savage fc Mathisl 



Il979l ) for the Johnson-Cousin passbands a nd Af435w = 
4.02 E{B - V), A F62 sw = 2.62 E{B - V") l|Sirianni et ail 



Marcolini et al. 



0) 



0.5 



0.0 



-0.5 



0.5 



0.0 



-0.5 



DdSph 



D90 



D30 



a) 




1.5 -1.0 
[Fe/H] 



-1.5 -1.0 
[Fe/H] 



2.0 -1.5 -1.0 
[Fe/H] 



Figure 2. Upper panels: abundance ratio [a/Fe] plotted against [Fe/H] of Ng = 1000 sampled stars for the three different regions as in 
Fig.l. Lower panels: mean values of [a/Fe] for the same stellar sampling. 



-0.5 



o 



N -1.5 



o 




-2.5 -2.0 -1.5 -1.0 -0.5 -2.5 -2.0 -1.5 -1.0 -0.5 -2.5 -2.0 -1.5 -1.0 -0.5 
[Fe/H] [Fe/H] [Fe/H] 

Figure 3. [Fe/H] of JVg = 1000 sampled stars plotted against the metallicity Z for the three different regions as in Fig. [T] The dashed 
line represents the linear relation given by eq. (2) in the case of pure SN II enrichment. 



120051 ) for the ACS passbands. For each star a random pho- 
tometric error extracted from a Gaussian distribution with 
a equal to the average photometric error estimated in that 
magnitude range by the above authors has been added. 



2.4 Adapting the model to ui Centauri 

The 3D simulations employed here were generated using the 
Marcolini et al. (2006) hydrodynamical code. As pointed out 
in Section [T] in the scenario that assumes ui Cen is the rem- 
nant of a dwarf galaxy accreted by the Milky Way, most of 
the initial mass of the proto-galaxy has been stripped, pre- 
sumably by the tidal interaction with the Milky Way. Un- 



fortunately, our model cannot account for the self-consistent 
interaction of the proto-galaxy with the Milky Way. Thus, 
we adopt a model ori ginally tailored to the Draco dSph (see 
iMarcolini et al.1 2006). but now focusing our analysis on the 
central regions of the models - i.e., those regions expected 
(potentially) to survive as a "globular cluster" like u) Cen. 

In the model, the SFH is represented as a sequence 
of instantaneous star bursts separated by quiescent periods 
(At = 60 Myr). At the end of the entire star-formation 
process, nearly 6 x 10 5 Mq of stars have been formed, and 
6 x 10 3 SNe II exploded. The Draco galaxy has a stellar 
mass content smaller by a factor of 4-8 (e.g. iMevlan et al.l 
Il995l ; lvan de Ven et alj|2006l ; iMateoll 19981) and a much more 
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Figure 4. Mctallicity (Z) distribution function of the long lived stars for the "entire dSph volume" (left panel), the "central region" 
(middle panel) and the "nuclear region" as in Fig.[T] The dashed lines represent stars with [Fe/H]J^ —1.4. 



extended stellar distribution than u Cen. These differences 
are possibly even more marked if we consider the structure 
of u> Cen before the interaction with the Milky Way. Thus, 
in order to match the chemical properties of the dominant 
population of the system, the number of SNe and the initial 
amount of gas in the "proto" uo Cen have been scaled appro- 
priately, and the star formation suppressed after ~1 Gyr. 

u) Cen has a stella r mass content of ~ 3 x 10 6 M© 
van de Ven et all ([2006) and a radial profile c onsistent with 
a Kin g model with r c = 4.1 pc and r t = 83 pc l|Ferraro et al.l 
2006). We followed the chemical evolution of the stellar pop- 
ulation of the dSph progenitor in three different regions,, 
defined as follow: i) the "entire dSph volume" with a tidal 
radius r*,t = 650 pc (see Marcolini et al. 2006 for more 
details); ii) the "central region" within a radius r*, t = 90 
pc; and, Hi) the "nuclear region", with r*^ = 30 pc. These 
latter two regions are dense enough to survive the tidal dis- 
ruption incurred through the interaction with the Galaxy, 
and appear today like a globular cluster. We thus focus our 
analysis on the central and nuclear regions, and will refer to 
them as models D90 and D30, respectively. 



3 RESULTS 

Figure [T] illustrates the [Fe/H] distribution of the long-lived 
stars (with masses M < 0.9 M©) formed over 1 Gyr in the 
three regions described above. The distribution shown in the 
left panel is a rather smooth log nor mal-type distribution , 



similar to those f ound i n dSphs (e.g .l Bellazzini et al.l 120021 
Babusiaux et al.| 20051; iKoch et alj [ 
L||2007i) 



2006: 



Battaglia et al.l 
The profile in the "rentral re- 
shows a bimodal structure 



20061 : iBosler et al.l |200" 
gion" (middle panel), though 
similar to that observed in w Cen (e.g. iNorris et al.lll996l : 
ISuntzeff fc KraftJll996h . with a maximum at [Fe/H] = -1.6 
and a secondary peak at [Fe/H]~ —1.3 accounting for ~ 25% 
of the cluster's stellar content. The different shape of the 
latter's distribution is due to the higher SN la rate in this 
region, as well as to the fact that external iron rich SN la 
pockets converge toward the centre during the re-collapsing 



stages of the ISM evolution. The right panel in Fig. \T\ shows 
the distribution in the "nuclear region". The noise present 
is due to sampling statistics and to the fact that the vol- 
ume considered is comparable to that of single pockets. In 
each panel only a negligible fraction (< 1%) of stars has a 
metallicity [Fe/H]> -0.7. 

The upper panels of Fig.[2]show the [q/Fc]-[Fc/H] plane 
for N s = 1000 sampled stars. As in Fig. [1] the left panel 
refers to the whole galaxy while the central and the right 
panels refer to the central and the nuclear regions, respec- 
tively. The lower panels represent instead the mean value 
of [a/Fe]. Because of the SN la contribution, the Fe-rich 
stars are, as expected, a-depleted: we find mean values of 
([a/Fe]) ~ 0.4 at [Fe/H]= -1.7 and {[a/Fe]) ~ 0.2 at 
[Fe/H]= —1.3, with a large spread in the [a/Fe] distribu- 
tion at [Fe/H]> -1.5 . 

In general, our mean value of [a/Fe] in the high metal- 
licity regi me (Fe/H> —1.0) i s slig htly smaller than th at ob- 
served bv lSmith et~al] (|2000i ) and lOriglia et al.l (120031 ). This 
discrepancy is most likely due to the simple form of the 
SFH we assumed. In Section [3TT] we will show that this dis- 
crepancy, as well as the lack of metal-rich ([Fe/H]> —0.7) 
stars, may be reduced by adopting a more protracted SFH. 
Indeed, assuming that the star formation can last beyond 
1 Gyr at a lower rate, high-metallicity stars ([Fe/H]> —1.3) 
may continue to form with high [a/Fe] values because of the 
freshly-expelled ejecta from recent SNe II explosions. 

We must re-iterate the difference between the pollution 
of iron and that of metals, as a whole. It is clear from Table 1 
that the ratio between the iron mass and the mass of all the 
metals is rather different in the ejecta of the two types of 
SNe. Thus, when the contribution of the SNe la to the ISM 
pollution becomes relevant, the metal mass fraction Z is not 
linearly proportional to the amount of iron [Fe/H]. It can be 
shown that, with our assumptions, it is: 

Z = 0.015 xl0 [Fe/H1 10 [Q/Fel [1 + 0.08 (2.82 x W~^ F ^ - l)] 

Figure [3] shows the log(Z/^Q)-[Fe/H] diagram for the 
same sample of stars of Fig. [2] The dashed line represents 
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Figure 5. Upper panel: temporal profile of tha assumed star for- 
mation rate in the LONG model. Lower panel: age-Z and age- 
[Fe/H] distributions of 1000 sampled stars for the D90LONG 
model. Note the metallicity spread among coeval stars. 



the log(Z/ZQ)-[Fe/H] relation for stars formed in regions 
enriched only by SNe 110: 

log(Z/Z ) = O.42 + [Fe/Jf], (2) 

assuming Z Q = 0.016 (|Grevesse fc Sauvallll998h . 

The Z distribution for the long-lived stars is plotted 
in Figure [4] The bimodal structure of the [Fe/H] distribu- 
tion (see Fig. [1} is now absent and the spread of the Z 
distribution is similar to the spread of the dominant [Fe/H] 
poor population. In fact, the number of SNe la occurring 
after 1 Gyr accounts for ~ 20% of the iron content but only 
~ 1.5% of the produced metal mass. Therefore, in the re- 
gions polluted inhomogeneously by SNe la, the [Fe/H] ratio 
increases dramatically, but the metal mass fraction Z re- 
mains basically unaffected (c.f. Figure This effect is par- 
ticularly important in the central regions, where the stars 
formed in "SNIa pockets" can constitute ~30% of the entire 
cluster population. In Figure [4] the Z distribution for stars 
with different iron content is shown for the three regions 
described in Section [2.41 Note that, in the central regions, 
Fe-rich ([Fe/H]> —1.4) and Fe-poor stars span nearly the 



1 This relation is obtained from the previous formula assuming 
[o/Fe]=+0.45, representative of SNe II ejecta (cf. Table 1) 



same Z range. As we will see in Section |4j this difference 
influences the position of the stars in the CMD. 

3.1 Longer star formation history 

As discussed in the previous section only a negligible frac- 
tion of stars in our model have a [Fe/H]> —0.7 and the mean 
value of [a/Fe] in the high metallicity regime ([Fe/H]> —1.0) 
is slightly smaller than that observed. These discrepancies 
can be traced to the simple form of the SFH adopted. In- 
deed, in Section \2. 41 we assumed that the star formation of 
our reference model is suddenly halted after ~1 Gyr, due 
to the interaction with the Milky Way. However, the ef- 
fect of ram pressure stripping and tidal interactions in more 
massive dwarf galaxies is likely to be more prolonged, and 
the time to completely strip the gas from such a system 
may be of the order of several Gyrs, depending upon the 
stre ngth of the ram pressure and the orbit of the system 
(e.g. iMarcolini et~aLll2003l ; iMaver et al]|2006h . 

To explore this possibility, we assumed that after 1 Gyr, 
the system still forms stars for a further 0.6 Gyr, but at a 
lower rate. This hypothesis is consistent with the ancient 
proto-galaxy losing most (but not all) of the gas during its 
first interaction with the Milky Way. This newly adopted 
SFH is plotted in the upper panel of Figure [5] In the follow- 
ing, we refer to this model as LONG. 

The reduction of the ISM reservoir during this 0.6 Gyr 
period leads to a corresponding reduction in the rate of the 
SNe II, whose progenitors are short-lived stars. Conversely, 
in the same period, the SN la rate does not decrease sig- 
nificantly because, in this case, the precursors have longer 
lifetimes and their explosion rate still depends on the past, 
higher, SFR. Thus, the rapid gas depletion, together with 
a rather unaffected iron production, leads to a rapid metal 
enrichment of the stars forming after the stripping episode, 
during the final 0.6 Gyr. 

As already discussed in Section 12.41 in our original 
model all the starbursts had the same intensity, and the 
ISM remains entirely bound to the galaxy. Unfortunately, 
we cannot mimic the gas stripping by simply reducing ar- 
tificially the amount of ISM in our model; in this way, in 
fact, the metals mixed with the ISM would be removed as 
well and the relative role of SNe II and SNe la as polluters 
would not change. Instead, the metal enrichment is obtained 
by increasing in time the quantity of SN ejecta present in the 
simulation; in order to take into account the corresponding 
drop of the SNII rate due to the decreasing SF over the last 
0.6 Gyr, in parallel with the rather unaffected SNIa rate (see 
above) , the relative rates of increasing metal enrichment due 
to the two SNe types are chosen a posteriori to best fit the 
data. Not having strong constraints on the orbit of first in- 
teraction of the proto-galaxy with the Milky Way (and how 
the gas stripping proceeded), this remains a necessary free 
paramenter. 

In the lower panel of Fig. [5] we plot the age-metallicity 
diagram for the D90LONG model. As evident, the metal- 
licity increases with time even if a large spread is always 
present, in agreement with the most recent determinations 
(e.g. Hilker et al. 2004; Stanford et al. 2006; Villanova et al. 
2007). Note that the spread is more evident in the [Fe/H] dis- 
tribution, as opposed to the Z one, because of the stronger 
influence of the inhomogeneous SNe la enrichment on the 
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Figure 6. Left panels: theoretical (D90LONG model) and observative i Sollima et_alj ]2005a ) [Fe/H] distributions of long lived star 
(upper and lower panels, respectively). Right panels: theoretical and observative jNorris et alTll996h [Ca/H] distributions of long lived 
stars (upper and lower panels, respectively). 



Fe abundance. We find that we can generate the metal- 
licity distribution functions (and associated observational 
constraints, as noted below) of ui Cen assuming a SFH last- 
ing 1.6 Gyr, i n agre eme nt with the timescale proposed by 
iFerraro et at] (120041 ) and lSollima et alj l|2005bl ). 

In Figure|6]we show the [Fe/H] and [Ca/H] distributions 
obtained by the D90LONG model, together with those ob- 
served. Comparing the [Fe/H] distribution in the upper-left 
panel with the corresponding one in Figure [1] we note that 
the shape of the [Fe/H] distribution remains unaffected at 
low metallicity ([Fe/H]< -1.2), but for [Fe/H]> -1.2 the 
D90LONG model predicts a larger number of stars and a 
small peak at [Fe/H]~ —0.6. The agreement w ith the [Fe/H] 
distribu tion inferred by the observations of ISollima et al.l 
(2005a) is now quite good. 



In the same figure we also compare the [Ca /H] distribu- 
tion fo und by our model and the one inferred bv lNorris et al.l 
(1996). Again, the two distributions are in good agreement. 
Actually, we tend to overestimate the number o f stars having 
[Ca/H]> -0.8, but as stated by iNorris et al.1 |l996l) . their 
selection criteria tends to underestimate the number of stars 
in this range. Note that, contrary to the iron, the amount of 
calcium ejected by a single SN la is comparable to that ex- 
pelled by a SN II (SN Ia ej , C a/SN II ej , C a ~ 1.4fl Thus, the 
two peaks visible in the [Fe/H] distribution become nearly 
indistinguishable in the [Ca/H] distribution. 

In Fig. [7] we plot the projected radial distribution of 



Note that this ratio is ~ 8 for iron 
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Figure 7. [Fe/H] (upper panel) and [Ca/H] (middle panel) pro- 
jected radial distributions of 500 sampled stars for the D90LONG 
model. The observed [Ca/H] projected radial distributions by 
Norris et al. (1997) is shown in the bottom panel. Radius is given 
in pc. 



[Fe/H] (upper panel) and [Ca/H] (central panel) for a sample 
of 500 stars. As already discussed, the most iron-rich stars 
are more concentrated toward the centre (c.f. Section [2,1[) , 
In the bottom pan el, the [Ca/H] distribution observed by 
iNorris et al.l (|1997m is shown. We find good agreement be- 
tween the observed distribution and that predicted by the 
D90LONG model. 

Figure [8] (left panel) illustrates the log(Z) distri- 
bution predicted by the D90LONG model. The peak at 
log(Z/ZQ) = —0.65 is due to stars born in the final episode 
of star formation which have —0.8 ^[Fe/H]$; —0.4. These 
stars have ratios [a/Fe] in the range and —0.3 ^ [a/Fe] ^ 



+0.3, as shown by the two central panels of Figure [8] Com- 
paring the central panel of this latter figure with the cen- 
tral panels of Figure [21 we note that the mean [a/Fe] val- 
ues at [Fe/H]=-1.7, -1.3, and -0.6 are now < [a/Fc]>~ 
+0.4, +0.2, and +0.0, r espectively, in good agreement with 
the values inferred by lOriglia et al] ([2003). The improve- 
ment at higher [Fe/H] is due to the temporal tail of the star 
formation. Indeed, as the star formation proceeds (though 
at a lower rate), the new SNe II explosions release fresh a- 
elements which raise the [a/Fe] ratio of most of the newly- 
formed stars. Finally, the right panel of Fig.[8]emphasises the 
dispersion of the stellar metal fraction in the log(Z)-[Fe/H] 
diagram. 

To better compare our model with observations we 
plot separately in Figure [9] the abundance patterns of 
several a-elements (O, Mg, Si, and Ca) as predicted 
by our model, together with a dataset collected from 
the l i terature ( Francois et al] 1 19881; Brown fc Wallersteinl 



19931; INorris fc Da Costal Il995l; ISmith et al] Il995l. 



2000 



Pancino et al.ll2002l; IVanture et al.ll2002l ; lOriglia et al.ll2003l : 

Villanova et al] 2007 ) . Altough the data are in good agree- 
ment, discrepancies as large as 0.5 dex, are occasionally 
present (for an ex tensive discussion of this dataset see 
iRomano et"ai]|2007j l 

For each a-element we are able to reproduce the gen- 
eral trend of a nearly constant value up to [Fe/H] ~ — 1.0, 
followed by a slight decrease down to the value [a/Fe]~+0 
at [Fe/H]~— 0.6. We also simulate satisfactorily the spread 
of the data and the peculiar stars with sub-solar [a/Fe] (note 
in particular the agreement between observations of [O/Fc] 
at [Fe/H]> —1.6 and the model prediction). 

An alternative explanation for the O and Mg deple- 
tion is the formation of stars from gas polluted by AGB 
ejecta (c.f. ICottrell fc Da Costalll98ll ). In this contest, self- 
consistent models of the c hemical evolution of the globular 
cluster NGC 6752 made bv lFenner et all l)2004h are not able 
to reproduce the O depletion, while Mg turns out to be pro- 
duced rath er than destroyed; see a lso lDenissenkov fc Herwid 
l|2003h and lKarakas et all (|2006bh . 

The uncertainties that affect the underlying AGB mod- 
els undermine the reliability of the predictions and may le ave 
room for an AGB solution l|Ventura fc D 'AntonalliooBal lbl). 
The choice of convection treatment within the AGB mod- 
els is particularly important. When convection is mod- 
eled more efficiently, deep Hot Bottom Burning can result 
ijBloecker fc Schoenbernerlll99ll . e.g.), with a strong deple- 
tion o f oxygen and magnesium (e.g. IVentura fc D'Antonal 
l2005bl la). while the C+N+O remains constant to within a 
factor of two. This could help to explain the presence of low 
oxygen abundance stars in the top right panel of Figure 10. 

Since intermediate mass stars (4-8 M e ) evolve over 
timescales (30-70 Myr) comparable to that of the peak of 
SNe la activity (in a "burst" environment) we expect that 
regions inhomogeneously polluted by SN la should also be 
enriched by intermediate-mass AGB stars, complicating the 
picture. Since AGB stars produces large q uantities of helium 
(up to Y=0.35, e.g. lD'Antona et al.ll2005r ) one might expect 
stars polluted by AGB ejecta to be rich in helium, even if 
the content hypothesised by s everal authors to explain the 



mult i ple populations in u > Cen ijPiotto et al.ll2005T: iLee et al 



2005: lSollima et al.l 



20051 ; IPiotto et al] 



2007) and in other GC (|D'Antona et al 
2007) require that such stars formed 
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Figure 8. Left panel: log(Z) stellar distribution. Central panels: [c*/Fe]-[Fe/H] diagram for 1000 sampled stars (top), and ([a/Fc])-[Fc/H] 
relation (bottom), where the brakets indicate the mean value of [a/Fe]. Right panel: log(Z/ZQ)-[Fe/H] diagram of 1000 stars sampled 
in the central region. All the results refer to the model D90LONG 



from "pure" AGB ejecta. Regardless. iKarakas et al.l (|2006al ) 
found that only a maximum helium content of Y~ 0.26 was 
feasible whitin their models, without violating the observa- 
tional constraint of near-constant C+N+O. 

Finally, lo Cen shows a clear trend of increasing heavy 
elements produced by s-process s ynthesis (e.g. Rb, Ba, La , 
and Nd) with inc reasing [Fe/H] (|Norris fc Da Costal fl995; 
ISmith et al.l l200Cjh . This is usually interpreted as a sign of 
the progressive chemical enrichment from 1.5-5 Mq AGB 
stars. The stellar lifetim e of stars with mas s in these range 
(~ 1 Gyr for a 2 M Q star lSchaller et al.lll992T ) are compatible 
with the SFH inferred by our study and furthermore the 
mass reduction in the last ~600 Myr can aid in the rapid 
enrichment of these long-lived progenitors' elements. This 
kind of s-process enrichment trend is peculiar to lo Cen, in 
comparison with other GC stars. 

To be fair, models of the sort described above may only 
apply to "traditional" GCs, and may not be applicable to 
lo Cen, but we felt it prudent to at least discussion the fea- 
sibility of this scenario here. Recall, lo Cen is the only GC 
clearly showing a significant spread in iron abundance (as 
required by SNe la enrichment) . We are currently developing 
a new model for canonical GC evolution based upon similar 
ingredients (Marcolini et al. in preparation). 

In conclusion, AGB pollution has certainly occurred at 
some level within lo Cen, and cannot be ignored. With this in 
mind, in the next section we test the validity of our model 
by simulating the CMD of lo Cen and comparing it with 
observed CMDs from the literature. 



4 COMPARING THE COLOR-MAGNITUDE 
DIAGRAMS 

Most of the evidence gathered in the past regarding the 
metallicity spread, the presence of multiple populations, and 
the possible helium enhancement in lu Cen have been de- 
duced from the peculiar morphology of different sequences 
in the CMD. It is interesting to check how our models com- 
pare with the most recent and accurate observational data. 
As pointed out in Section [3] our model predicts a wide 
spread in the distribution of [Fe/H] as well as [a/Fe]. It has 
been shown that both [Fe/H] and [a/Fe] produce similar ef- 



fects on the stellar isochrones (ISalaris et al. 19931; fKim et al.l 
120021 ; ICassisi et al] | |2004 iPietrinferni' et al. 2006, although, 
see iDotter et al.l 20071 for a more thorough discussion of the 
impact of individual elements on stellar population models). 
In particular, the low first-ionization potential of Mg, Si and 
Fe favours the formation of H~ ions that are the main con- 
tributors to the mean opacity. Moreover, the increase in the 
abundances of C, N, O and Ne produces significant change in 
the opacity and in the efficiency of the CNO nuclear burni ng 
that mimic a higher average metallicity Z <|Renzinilll977h . 

In Fig. [10] the color distribution of RGB stars at two 
different magnitude levels is shown fo r both the simulated 
(D30LONG model) and the observed (ISollima et al.ll2005aT ) 
CMD in the V,B-V plane. The colour distribution of the 
RGB of lu Cen is well reproduced. The most metal-rich 
anomalous component (RGB-a, [Fe/H]~-0.6 ), well distin- 
guishab le in the red part of the CMD of ISollima et al.l 
(2005a|), is less defined in our model (although still present) 
because of the large spread in [a/Fe] ratios predicted at 
this metallicity. However, given the relative simplicity of the 
model, we consider the agreement between the simulated 
and the observed CMD satisfactory. 

A behaviour similar to that observed in Fig. [10] is also 
present in Fig. 1111 where our D30LONG model is com- 
pared with the ACS CMD by iFerraro et all (|2004T ) in the 
sub-giant branch (SGB) region. As can be noted, the ob- 
served magnitude spread of the upper part of the SGB, as 
well as the presence of the anomalous SGB (SGB-a), are ac- 
counted for. However, the shape and the slope of the SGB-a 
are not w ell reproduced . This effec t has been already di s- 
cussed by IFerraro et all (|2004T ) and ISollima et all i|2005bl . ; 
and references therein) and could be linked to the possible 
presence of a significant He overabundance at these metal- 
licities. According to our model, stars lying along this se- 
quence should have a high iron content ([Fe/H]~ —0.6), in 
agreement with spec troscopic measurements of stars belong- 
ing to this sequence (fkilker et al. l l2004l ; ISollima et alll2005bl ; 
IStanford et"al]|2006h . It is worth recalling that none of our 
previous models (DdSph, D90 and D30) reaches such a high 
metal content. Indeed, this anomalous metal- rich population 
appears to form in the last ~0.6 Gyr of evolution of lo Cen. 

Finally, in Fig. 1 121 we c ompa re our model with the CMD 
obtained bv lSollima et al.l l|2007h from deep FORS1 observa- 
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Figure 9. Predicted abundance ratios of several a-elements to iron as a function of [Fe/H] of 7Vg = 500 sampled stars (left panels). 
Observative data are in the right panels (see text for the references). 



tions of the MS of u Cen. Metal-rich stars ([Fe/H]> -1.4) 
are marked as grey points in the synthetic CMD. In the 
bottom panels, the colour distributions of Fe-poor ([Fe/H]< 
— 1.4) and Fe-rich ([Fe/H]> —1.4) stars in the simulated 
CMD and of observed stars are shown. As can be noted, Fe- 
rich stars are mostly located in the same region of the other 
Fe-poor stars as expected for two population with similar 
distributions of Z (see Section [3j . 

As noted in Section 1, spectroscopic analyses indicate 
that metal-rich ([Fe/H]~ —1.3) MS stars in uj Cen have been 
found to lie in the blue portion of the MS, in contradictio n 
with the canonical theoretical models dPiotto et alj |2005). 
In order to explain the anomalous location of these stars in 
the CMD, a large helium overabundance has been claimed 



by several authors ( 


Norrisjl2004lPiotto et al.ll2005l;lLee et all 


2005l;ISollima et al. 


2007). According to our model, the mean 



helium abundance remains basically unchanged during the 
entire process of self-enrichment at any distance from the 
cluster centre (AY < 0.05). The same conclusion is drawn 
from recent analytical (ie. without associated dynamical 
evolution of the ga s) galactic chemical evolution models 
|Romano et al.ll2007h . Although our model fails to reproduce 
the MS morphology of u Cen, the lower mean a-element 
abundance predicted for the Fe-rich group of stars tends to 
reduce the color distance bet ween bMS and rMS. A ccording 
to the theoretical models of IStraniero et all (|l997t ). in this 
scenario the required amount of helium overabundance to 
shift the metal-rich ([Fe/H]~ -1.3) MS on the blue side of 
the dominant Fe-poor cluster MS is reduced by a factor of 
40% (AF ~ 0.09). 

Finally, in Fig.QUwe plot our full (R,B-R) CMD. The 
location in the CMD of stars in three different ranges of 
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Figure 10. Comparison between the simulated (D30LONG model; left panel) and the observed (from Sollima et al. 2005a; central panel) 
(V; B-V) CMD of lo Cen. Right panels shows the color histograms between two different V magnitude intervals (indicated by rectangles 
in the left and central panels) for the observed and simulated (grey histograms) CMD. 



[Fe/H] are shown. Note that, although the three groups pop- 
ulate preferentially different regions of the CMD, they par- 
tially overlap. As discussed in Section [3] this is due to the 
fact that the [Fe/H] is not stictly linearly proportional to the 
metal mass fraction, and stars with the same Z can have 
different [Fe/H] content. In this picture our model could 
explain the peculiar position of the outliers observed by 
IVillanova etaD (|2007h in their spectroscopic analysis of sub- 
giant stars. These stars have been indeed found to have a 
[Fe/H]>-1.0 although their position in the CMD would sug- 
gest a significant smaller [Fe/H]. The comparison with our 



model suggests that the location of these stars in the CMD 
could be explained if these stars show a significant [a/Fe] 
depletion. 

While a detailed comparison of every aspect of the CMD 
is beyond the scope of the paper, we have demonstrated 
that the possibility that uj Cen had its origin in a single 
progression of star formation and metal enrichment cannot 
yet be entirely dismissed (cf. Villanova et al. 2007). Our 
proposed scenario is not meant to be interpreted as the only 
solution to the problem and to recover all the aspects of the 
CMD morphology, but neither should it be ignored when 
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Figure 11. Comparison between the simulated (D30LONG 
model; left panel) and observed (from Ferraro et al. 2004; right 
panel) (F625W; F435W-F625W) CMD of u Cen. 



Figure 13. Simulated (R; B-R) CMD of u Cen (D30LONG 
model). In the different panels stars with [Fe/H]<-1.4 (left panel), 
-1.4 <[Fe/H]<-0.8 (central panel) and [Fe/H]>-0.8 (right panel) 
are marked as black points. 



interpreting the MS, SGB and RGB morphologies (it may be 
that a required solution involving our mechanism alongside 
that of enhanced helium may be required). 



5 DISCUSSION AND CONCLUSIONS 

We studied the chemical evolution of the peculiar stellar 
system uj Cen under the assumption that it represents the 
nuclear region of a dSph deprived of most of its mass by 
its interaction with the Milky Way. We based our model on 
hydrodynamical simulations by Marcolini et al. (2006) de- 
scribing the evolution of an isolated dSph galaxy. We focused 
our analysis on the central region of the system, comparing 
its chemical properties with the observational data available 
in literature. 

We paid particular attention to the different role played 
by SNe II and SNe la in the chemical enrichment of the 
gas. While the SNe II pollute the ISM rather uniformly, the 
SN la ejecta may remain temporarily confined inside rela- 
tively small pockets of gas before being mixed with the ISM. 
The stars forming in such pockets have lower [a/Fe] ratios 
than the stars forming elsewhere which show a significant 
a-enhancement (([a/Fe]) ^0.35). a-depleted stars represent 
only ~7% of the entire dSph stellar population, but their 
fraction rises up to ~30% in the central region. This differ- 
ence depends on the hydrodynamical evolution, as well as 
on the higher rate of SNe la in the central regions of the sys- 
tem. As a consequence, the chemical properties of the stars 
forming in the nuclear region of a dSph differ substantially 
from those of the main dSph population, and can resemble 
those observed in u Cen. 

The Z and [Fe/H] distributions of u) Cen have been 
reproduced assuming that 85% of the stars were born within 
1 Gyr (in a sequence of equal-intensity starbursts) and the 
remaining fraction formed during the subsequent 0.6 Gyr 
at a lower rate (owing to the interaction with the Galaxy). 



The stars created in this last period show a rapid increase 
of their [Fe/H] because of the rapid reduction of the gas 
mass content as a consequence of the interaction between 
the dSph and the Milky Way. 

The first constraint against wich our model has been 
tested is that of the [Fe/H ] and [Ca/H] distr i butions of the 
stars observed in uj Cen l|Norris et al 1 ll996l : ISollima etajj 
l2005al ). The inhomogeneous mixing of the metals ejected by 
SNe la naturally explains the observed multimodal distribu- 
tions, as well as their radial profile. 

From accurate CMD simulations, we compared the mor- 
phology of the evolutionary sequences with that predicted 
by our model. Our model is successful in reproducing the 
overall spread in the RGB and SGB morphology, but fails 
to mimic the SGB-a morphology. It also fails to reproduce 
the bimodal MS observed in this cluster. However, due to 
the relationship between the [a/Fe] and [Fe/H] ratios in the 
SNe la pockets, the discrepancy between the predicted and 
observed location of the bMS component can be ameliorated 
by requiring a reduced level to the helium o v erabundance 
(A Y ~ 0.09) usually i nvoked l|Norrisl |2004| ; IPiotto et all 
l2005l ; ISollima et al.ll2007l ). 

A further test able to validate the inhomogeneous pol- 
lution expected by our model should be the observation in 
uj Cen of some stars with sub-solar [a/Fe] ratios (at [Fe/H]^ 
— 1.5). Note that stars with such a peculiar abundance pat- 
terns have already be en observed in the centre of some 
dSph s like Sagittarius ( Bonifacio et al. |2004|; [Monaco et al 



2005 ; iMcWilliam fc Smecker-Hane 120051 ; ISbordone etal. 
I2007T). Draco dShetrone et all 120011) F ornax and Sculptor 



( Venn fc Hil]||2005l ; iTolstov et al.ll2006l ). 

We note that in order to rea lise our synthetic 
CMDs, we utilised evol utionary tracks (|Cassisi et al.ll2004l ; 
iPietrinferni et a l. 2006) calculated for two different a- 
enhancement levels ([a/Fe]=+0.0 and [a/Fe]=+0.4). We 
then extrapolated down to the ratio [a/Fe]=— 0.4 (for each 




Z value) in order to take into account the entire [a/Fe] range 
covered by the stars in our model. Unfortunatelty, most of 
these stars do not have solar-scaled abundances, as assumed 
in the utilised isochrones. Tracks taking into account the 
relation among Z, [Fe/H] and [a/Fe] given by equation (1) 
would be desirable to obtain more realistic tests. Moreover, 
further models taking into account the dynamical interac- 
tion of the parent dwarf galaxy with the Milky Way and the 
effect of AGB pollution will be required. 
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